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Abstract. In the search for vacuum solutions, with or without a cosmological 
constant, A, of the Einstein field equations of Petrov type N with twisting principal 
null directions, the CR structures to describe the parameter space for a congruence 
of such null vectors provide a very useful tool. A work of Hill, Lewandowski and 
Nurowski has given a good foundation for this, reducing the field equations to a set 
of differential equations for two functions, one real, one complex, of three variables. 
Under the assumption of the existence of one Killing vector, the (infinite-dimensional) 
classical symmetries of those equations are determined and group-invariant solutions 
are considered. This results in a single ODE of the third order which may easily be 
reduced to one of the second order. A one-parameter class of power series solutions, 
g{w), of this second-order equation is realized, holomorphic in a neighborhood of the 
origin and behaving asymptotically as a simple quadratic function plus lower-order 
terms for large values of w, which constitutes new solutions of the twisting type N 
problem. The solution found by Leroy, and also by Nurowski, is shown to be a special 
case in this class. Cartan's method for determining equivalence of CR manifolds is 
used to show that this class is indeed much more general. 

In addition, for a special choice of a parameter, this ODE may be integrated once, 
to provide a first-order Abel equation. It can also determine new solutions to the field 
equations although no general solution has yet been found for it. 
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1. Introduction 

The search for general classes of solutions of the Einstein equations that either have pure 
vacuum for their source or a non-zero cosmological constant, that are of Petrov type N, 
and have principal null rays with non-zero twist has been continuing for a very long time 
now. We will use the (relatively common) nomenclature "Einstein spaces" for solutions 
with this sort of a source, i.e., either the pure vacuum or that vacuum with a non-zero 
cosmological constant. A, appended to it. With this idea firmly in mind, there are only 
two solutions known so far for twisting type N Einstein spaces: the one described by I. 
Hauser [HE], which has A = 0, and the one by J. Leroy [3], with a non-zero value for A, 
but which in the limit as A — )■ degenerates to a fiat solution, rather than one of type 
N. Because of the apparent difficulty of the problem, many different approaches have 
been used to attempt the finding of such a solution. With a requirement of one or more 
Killing vectors, the problem can be reduced to the solution of a single, nonlinear ODE, 
which has been produced in several forms by different authors [H El El [7] ; nonetheless 
this approach has produced no new solutions. Looking at the problem as a reduction 
from complex-valued manifolds via Plebahski's hyperheavenly equation [8] has produced 
no new solutions [9]. Therefore we were quite interested when we became aware of a 
different approach via a recent paper by Pawel Nurowski [11] , looking for exact solutions 
of this type with non-zero cosmological constant. 

Many of Nurowski's research articles use the fact that one can productively study 
(4-dimensional) Lorentz geometries which admit a shearfree, geodesic null congruence 
of curves by viewing the 3-parameter space that picks out any particular curve in the 
congruence, as a (3-dimensional) CR structure [12]. In [10], he and his collaborators use 
the CR function in such a structure to create a very appropriate choice of coordinates 
for a twisting type N Einstein space, and reduce the Einstein equations to a set of 
nonlinear PDE's for a couple of functions of three variables. Then in [TT], he makes 
a clever ansatz depending only on a single variable and discovers a particular twisting 
solution; unfortunately that solution turns out to be the same as the one mentioned 
above as first found by Leroy, as he notes in a more recent paper [13]. However, we were 
quite intrigued by the approach and have made some efforts to follow it through with 
the hopes of obtaining more general solutions of the equations in Nurowski's article. 

A spacetime of type N allows one, and only one, congruence of twisting, shearfree, 
null geodesies, referred to as a principal null direction (PND); this (3-parameter) family 
of null geodesies allows the option to choose a single coordinate r along any such geodesic, 
and to associate the other three degrees of freedom in the parameter space as a model 
for a CR manifold. We first insist that our manifold admit a Killing vector in the 
real direction in this (3-dimensional) CR manifold, so that the remaining unknown 
functions depend only on the complex coordinates there, and then calculate the (infinite- 
dimensional) classical symmetries for the system. This allows us to derive a quite simple 
nonlinear third-order ODE which the invariant solutions of the classical symmetries must 
satisfy. Because this equation does not contain the independent variable explicitly, it 
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can be immediately reduced to the following second-order ODE, for g = g{w), with two 
slightly different forms that differ by a constant: 

te:^_2c_^ ^^^^^^ 

We are then able to show that the Leroy-Nurowski solution is indeed a special solution 
for this equation. At this point it is worthwhile to introduce a question as to how 
one knows that the new solution of Nurowski does indeed describe localljll] the same 
manifold as the solution found by Leroy. The method was originally created by E. 
Cartan [151 HSl HI] to prove equivalencies of CR structures, without the need of actually 
determining an explicit coordinate transformation between the two sets of coordinates 
on two manifolds. Instead, one determines the values of a set of invariant quantities for 
a CR structure, the same for all equivalent such structures. Therefore it is necessary 
to calculate the invariants for Leroy's solution and compare them with the ones already 
known to Nurowski for his solution, noting that they are the same constants. We have 
therefore also calculated these invariants for our class of solutions, which we find to be 
quite different. 

For the case when the constant C in ([T]) is zero, an integral transformation may 
be performed to reduce that equation further, to a first-order ODE of Abel type, for 

/ = 



t \ 2J \ 3 J t \ 5/ 2t' 
which, quite unfortunately, we have not been able to identify as any of the known solvable 
types of Abel equations [TSl [H] . Nonetheless, we believe that these two equations, at the 
time of publication, are the simplest ODEs available that determine nontrivial twisting 
type N Einstein spaces. Returning to the case when C = 1, our current examinations 
suggest the possibility that the solutions of this equation might define a new class of 
transcendental functions, which constitutes the major result of this paper, establishing a 
new set of solutions to the type N problem with A ^ 0. The remainder of the paper will 
describe the process involved in this and our reasons for stating that these are indeed 
new solutions. In particular, we will present solutions to ([T]) in the forms of power series 
and Puiseux series, both shown to be locally convergent. Although prior to this time 
there were indeed only two known twisting type N Einstein spaces, it is always worth 
remembering that the solution space for the problem is in fact quite large. It has been 
shown by Sommers [20] that the full set of solutions for type N is given by two complex 
functions of two real variables. Surely the requirement of non-zero twist puts a very 
strong constraint on this, but it is expected that there should be a large number of new 
analytic functions involved in the full solution of the twisting type N problem. 



^AU our considerations are local, both in the Lorentz-signature spacetimc and in the associated 
complex spaces we need to use. 
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2. CR structures and reduced Einstein equations 

As already mentioned those (4-dimensional) Lorentz geometries which admit a shearfree, 
geodesic null congruence of curves have some history of people using CR structures to 
study them. They were first introduced into mathematics by Poincare and extensively 
studied by E. Cartan [151 HSj • Good sources of background on the relation of these two 
geometric concepts may be found, for instance in the thesis of Nurowski [21], and also 
in the very detailed discussion of their use for Einstein spaces in his article with Hill and 
Lewandowski [10], which generalized earlier work to include the possibility of a non-zero 
cosmological constant. 

The approach to a Lorentz-signature manifold begins with the usual form for the 
metric [22] in terms of a basis of 4 null 1-forms, modified for our choice of signature, 
(+,+,+,-): 

g = 2 {e^e^ + e^e^) , (2) 

where 6'^ is the complex conjugate of 6^ while 6'^ and 6^ are real, and the product is the 
usual symmetric product of 1-forms to create the metric. For Einstein spaces of type N, 
which have only one PND, we call its tangent vector k, which is a well-defined vector 
field in the neighborhood of each point, and take 6^ = g{k, ■) as its dual 1-form. The 
requirement that the twist of k should be non-zero is equivalent to the restriction that 
6^ A d^'^ ^ 0. We also choose a real coordinate r as a parameter along this congruence, 
and a real function P ^ on the manifold so that k = P~^dr {dr = d/dr, r G M). This 
congruence, still locally, is a 3-parameter family of null shearfree geodesies, so that at 
each value of r, the manifold still needs three additional coordinates. Following the work 
described in detail in [10], one is able to describe those degrees of freedom as coming 
from a 3- dimensional, strictly pseudoconvex CR manifold, M, equipped with (at least) 
one non-constant CR function (, and classes of pairs of 1-forms A (real) and /i (complex) 
such that 

A A /i A /i 7^ 0, /i = dC, = dC, 

dX = ifi A Ji + (c/i -|- c/i) A A 

with c a complex-valued function on M. A CR structure is a 3-dimensional real manifold 
M equipped with an equivalence class of pairs of 1-forms, (A, /i), as above. Another pair 
(A',/i') is considered equivalent to (A,/i), and therefore simply another representative of 
the same equivalence class, iff there are functions / 7^ (real) and h 0, g (complex) 
on M such that 

A' = /A, fi' = hfi + gX, ]j! = hji + gX. (4) 
This CR manifold can be lifted into a spacetime, so that the tetrad may be displayed 



in the following way [TO] : 



(5) 

= px, e^ = p (dr + Wfi + Wfi + HX) . 
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Of course this is exactly the form that a spacetime with such a distinguished PND is 
supposed to have [22] . 

Referring back to ([5]), we can see that the twist is just proportional to A A dA. 
Using (|3]) to determine this, we see that these forms have been chosen to put the non- 
zero character of the twist very clearly in evidence as non-zero. However, still looking 
at those equations, the function c that appears there is quite an important new function 
on the manifold. Using the closure of the equations in ([3]), one determines an important 
reality condition on the derivatives of c: 

dc = Be, 

once a dual basis of vector fields is introduced, which, however, is not a commutative 
basis: 

(9o,9, 9) dual to (A,/i, /i), 

[d,d] = -i(9o, [do,d] = cdo, [do,d] = cOq. 

At this point one has sufficient information to write out the explicit forms of the Einstein 
equations, which, in generality, say that R12 = A = with all the other components 
of the Ricci tensor vanishing. We quote from [10] and [H] which show that the results 
are the following: 

cos(|) 

where the functions a, q (complex) and /i, p (real), all independent of r, satisfy 
a = c + 291ogp, 

2 2 Bp — p idBp + Bdp) i 
q = -Ap + — ^0 logp - (9c, 

,2 



h = 2Ap^ + ^9pBp-p{dBp + Bdp) _ 2 

p2 



Given all the above, the functions a, c, h,p and q define a twisting type N Einstein space, 
of the form given in ([2]) and ([5j), iff the unknown functions c and p satisfy the following 
system of PDEs on M 

dc = Be (6) 

[dB + Bd + cd + cB+^cc+^ {dc + Bc)]p= |Ap^ (7) 

^3 = 0, i?33 = 0. (8) 

as well as one inequality 

$47^0, 

in order that the spacetime should not be conformally flat. In terms of those variables 
already deflned, the Ricci tensor component R33 and the Weyl scalars \E'3 and take 
on quite nasty-looking expressions 



R 



33 



— ((9 + 2c) {p'^dl) - 8A(^|Ap2 + 6 (c(9 + cd) logp 
-|- 129 log p Blogp + 3cc — Be — 2ido logpj cos"^ ^- 
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^3 



2i 

p2 



a/-4iA (291ogp + c) 



e''-/2cos3 



(d + c) {2d\ogp + c) +2{2d\ogp + cy 



6-'"/^ cos^ 



p^ 3 
where the function / is defined by 

I = d {dlogp + c) + {dlogp + c)^ . 

Following the procedure of [10] (with the use of Maple) which calculated the curvature 
tensor using Cartan's structure equations, we present the calculated \l/4 above {simplified 
with the use of \&3 = 0) for A 7^ 0, which we believe has not been previously published. 
Despite the frightening appearance of -R33, the equations ([7]) and \E'3 = together do 
imply the requirement -R33 = 0. This tells us that within the established formalism the 
twisting type N solutions to the Einstein equations automatically satisfy the condition 
for an Einstein space, i.e., vacuum with or without a cosmological constant. For A = 0, 
the statement is obviously true (see also [22] p. 451) and was used in [10] to prove the CR 
embeddability of twisting type N vacuums, without cosmological constant. For A 7^ 0, 
one uses ([7]) to substitute the term |Ap^ in R33 and notices that the resulting expression 
is a linear combination of d'^s and "^3. The equation i?33 = is therefore superfiuous 
for the type N problem, which facilitates our calculation greatly. 



3. Killing vector in the n-direction 

It is useful to understand the meaning of the operator d by introducing a real coordinate 
system {x, y, u) on M such that we have 

C = X + iy, dc_ = \{d^- idy) , ^ du + LdC + I^dC . . 

d = d^- Ldu, do = i{dL - dL)du, ~ i{dL - dL) ' ^ ' 

with L = L{(,(,u) a complex-valued function [23|. In addition, the function L relates 
to the function c in the following way 

c= -d\n{dL-dL) -d^L. (10) 

One major difficulty of fully solving the system (j6]|H]) is that unlike ordinary 
coordinate differentiations, the selected (dual) basis for the tangent space is not 
commutative, and, even worse, the operator d itself involves the unknown function 
c. When this dependence on the real coordinates is written out explicitly, the original 
PDEs become formidably lengthy. Instead of facing this entire conundrum, we have 
decided to circumvent it, at least in this paper, by looking at the special case where the 
unknowns p and c have no M-dependence, i.e., dop = = doC. Geometrically speaking, we 
insist that the spacetime admits a Killing vector in the u-direction. Such an assumption 
simplifies the problem greatly in that one can treat the operator d the same as d^, when 
acting on either p or c. This is a generalization of the assumption made by Nurowski 
[TT] . where it was simplified to just dependence on y, i.e, two Killing vectors assumed. 
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Theorem 1 {CR emheddahility [H]) A CR structure <^ with c = c(C, C) is CR 
embeddable. From this particular form of c, a u-independent form of the function L 
can be constructed as 



^(C,C) = -i j «(C,C)dC (11) 

with a real- valued function a 7^ satisfying 

dc^a = —ca, d^a = —ca. (12) 
Associated to this L, the tangential CR equation drj = yields a second CR function: 

r^ = u + i [[ a(C,C)dCdC. (13) 



Proof: Because of the restraint c?^c = d^c, the system ( fT2l) is compatible and has 
a real solution a 7^ 0. Henceforward one can directly check that ( ITT]) satisfies ( ITOj) and 
that f lT3l) satisfies the equation 9?7 = (9^ — Ldu)'r] = 0. Clearly, the CR functions rj and 
C are functionally independent, i.e., d( A drj ^ 0. Therefore we acquire a second CR 
function. 

For a given function c = c(C, C); the equation ( fTOi) . viewed as a PDE for L, may give 
rise to multiple choices of the function L, hence various A's. However, such an ambiguity 
only constitutes different representatives of the same CR structure. To see this, one may 
look into the six Cartan invariants (details in Section [6]) and notice that they are all 
uniquely determined by the function c = c(C, C) (see, e.g., fl28l) and Appendix B), given 
that the function r defined in ( l28l) does not vanish. For CR structures with r = 0, 
they are all locally equivalent to a 3-dimensional hyperquadric inside [151 113 • 
alternative proof would be to show that there always exists a coordinate transformation 
u — u{(X,u) {du/du ^ 0) that takes a function L = L{(X,u) satisfying ffTOj) to 
the M- independent form ( fTTl) . This can be confirmed by checking the compatibility of 
PDEs regarding such an existence. In conclusion, the CR structure on M is uniquely 
determined once a function c = c(C, C) is given. 

The converse of the last statement above is, however, not true. In fact various 
choices of the function c may correspond to the same CR structure. We will see examples 
of this in later sections. Related to this matter, our assumption of the function c being 
M-independent is thus not a CR invariant property. A function c = c(C, C) ^lay acquire 
M-dependence through the transformation (jl]) that takes one representative (A, fi) of the 
CR structure to another. 

Now we apply the assumption and the following notations 

df d(^f = fi, Bf d-^f = /2, / = /o, (/ = P, c and c only) 

and then rewrite the system (EHHD as 

Cl = C2, (14) 
2pi2 + CoPl + C0P2 + ^CqCoPo + i(ci + C2)po = |Apo, (15) 

P0P122 - P1P22 + 2coPoPi2 - 2copiP2 + 2cipoP2 + (ci2 + 2coCi)po = 2A(2p2 + coPq)pI, (16) 
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where the last equation arises from \l/3 = 0. These are the PDEs we aim to solve. 
Moreover, the Weyl scalar \E'4 reads 

4 e~'^/^ /r\ 

= -A [2poP22 + Qpl + IOC0P0P2 + (C2 + ^cI)pI] cos=^ j . (17) 

4. Infinite-dimensional classical symmetries 

We follow the standard procedure (e.g., [21]) to calculate the classical symmetries of 
the system fll4til6p . Since f|T6l) is generally complex, we have to include its complex 
conjugate as well in the calculation. Moreover, we treat ( fT4|) as a constraint and encode 
it and its differential consequences directly into the choice of intrinsic coordinates so 
that this equation no longer needs further attention. This gives us three PDEs for three 
dependent variables p, c and c which depend on two independent variables C and C- The 
intrinsic coordinates within the first four jets that are relevant to the calculation are 
chosen as follows 

PO) Co, Co, 
Pl,P2,Ci,Ci,C2, 
Pll,P22,Cll,C22, 
Plll,P222,Clll,C222- 

The rest of the jet variables, such as pi2, P122, Ci2(= C22) etc., can be expressed in terms 
of the intrinsic coordinates through the PDEs and their differential consequences. 

With a considerable amount of manual work on the algebraic computer program 
Maple, we have managed to find the classical symmetries with the generating section 
given by 

^ = - \{d(^A + d-^A)pQ - Api - Ap2, 
e = dlA- {d(;A)co - Aci - Aci, 
e = dfA - {d^A)co - Aci - Ac2, 

where A = A{() is an arbitrary function of ( that is sufficiently differentiable. The Lie 
bracket of two symmetries with, respectively, Ai{() and ^2(0 yields a third symmetry 
with a new A^i^Q given by 

As = [Ai,A2] := Aid<;A2 - A2d^Ai. 

Therefore, we indeed obtain an infinite-dimensional set of classical symmetries for the 
system ( IT^llT6l) . In particular, they reduce to translational symmetries for nonzero 
constant A, and scaling symmetries for A oc (. 
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5. Group-invariant solutions and reductions to ODEs 

Setting the generating section (\[^, 6, 6) to zero and solving for p, c and c, we are able 
to obtain a remarkable ansatz for the Einstein equations: 

,(C,c)^^, .,c,c)^M±M£)±^, fciffll^ (18) 

with a new real argument 

idC-/idc)=Im/|dC. (19) 

where the constant Ci and the undetermined functions Fi 2(2;) are all real- valued. One 
may easily verify these expressions by direct calculation. 

Substituting the ansatz into (fT5l) and (fT6l) and noticing that all dependence on 
A,A^O, except those in the argument z, can be factored out, we have a neat reduction 
from the PDEs to a system of two ODEs for Fi and F2 only: 

= -F[' + F2F[ + lAF^ - - 3F^ + Cl)Fi, 
= -F['Fi + {F[f + AF^ + F^F^. 

The satisfaction of the second ODE above is given by the introduction of a single new, 
real-valued function J = J{z) such that 

Fi = ±^^I, ^2 = ^ - AJ. (20) 
Then the first ODE simply becomes 

J'" = _ 2Ajr - — A(J')' - 2(AV2 + Cl)J'. (21) 

2(7 ' 3 

Since this ODE does not have the argument z appearing explicitly, we can lower the 
order of the ODE through the standard transformation 

/ = p(j) =^ J" = pp' =^ J'" = p{ppy 

and obtain an even simpler equation of the second-order 
_ _ (P- + 2AJ)' _2Cl _W 

2P P 3 ^ ' 

A solution P = P{J) to (|22I) can give rise to a solution J = J{z) to (12T!) at least locally 
by inverting 

z + Co = J -pi^d J (23) 

with Co constant. This solution will be physical if it also makes Fi 2{z) real- valued via 
(|2U]) . which requires that locally 

P(J) > 0, J' > and J real-valued. (24) 

Therefore we are only interested in solutions for J[z) that are monotonically increasing, 
or equivalently, positive P{J). 
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We can also consider the special case of ([22]) with Ci = and A 7^ 0, i.e., 

2P 3 ^ ^ 

By introducing the following integral transformation 

J = i exp Q fit) dt^ , P( J) = ^ exp (^2 J f{t) dt^ , 

of which the inverse has the form 

P , AJ2 



AJ2' JP'-2P' 
we can further reduce fl25|) to an Abel ODE of the first kind [28], as already noted in 
the Introduction: 

Once the general solution / = /(t,C2) is acquired with a constant C2, we can find the 
general solution P{J) of (125!) by solving the following ODE 



>J2' V JP'-2P' 
of which the solution is given by 

rZ/A 

P{J) = Z{J)J\ with = -lnJ+ / /(t,C2) dt + Ca. (27) 

Simple as both f[22|) and f[26|) may appear, so far we have had no luck finding their 
explicit general solutions. For more comments on ( 126|) and Abel ODEs in general, see 
Appendix A. 



6. CR equivalency as classical symmetry 

To identify new twisting type N Einstein spaces obtained from f[2T|) . we refer to the 
following theorem as a natural way to classify metrics equipped with CR structures. 

Theorem 2 ([10], see Theorem 1.2 and references therein) Let (A^,(yf) be a 4- 
dimensional manifold equipped with a Lorentzian metric and foliated by a 3-parameter 
family of shearfree null geodesies. Then Ai is locally a Cartesian product = M x M. 
The CR structure (M, (A, /i)) on M is uniquely determined by [Ai^g) and the shearfree 
null congruence on Ai. 

By definition, a type N spacetime at each point has a unique PND. In the case of 
vacuums (with or without A), a PND must be geodesic and shearfree [25]. Thus for 
every twisting type N Einstein space, the shearfree null congruence is unique. Hence 
to confirm a new twisting type N vacuum metric, it is sufficient to show that its CR 
structure is distinct from the one of known metrics. This can be routinely done by 
computing the six Cartan invariants [151 UHl E], which are denoted respectively by 

ai,9i,rij (complex), 

P 1,11X1 (real). 
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Cartan showed that two local CR structures are equivalent iff their six CR invariants 
(defined when r 7^ in fl28l) ) are identical, except possibly for a sign difference in both aj 
and rji [15]. With the assumption of w-independence, we can write down, for instance, 
the simplest invariant computed from the 1-forms defined in ([3]): 



«/(C,C) 



5rd(^r + rd(^r + 8crr 



r 

' 6 



V^(rr)7 ' (28) 

I [d^l + 2cl) , / = —d(^d^c — cd^c. 

Here the function r 7^ 0, following the notation of Cartan, is not to be confused with the 
coordinate r along the null congruence. For our calculated /9/,7/ and 9i, see Appendix 
B. Note that this «/ only relies on c{(, (), c and their derivatives, which is also the case 
for all the other Cartan invariants. We first point out a remarkable feature of these 
invariants computed from the ansatz (I18|19l) . 

Proposition Given the ansatz f fTSfTOj) . all the following quantities are independent 
of A{Q and A{Q except those in the argument z: 

In another word, they are all functions of z only, e.g., C) = Pii.^)- 

Proof: Except for a lengthy but straightforward symbolic computation with Maple, 
we are, at the moment, still not aware of any other more insightful way of proving this 
result. Here we only emphasize that the law ^Jv^Jw = ^Jvw is in general not true in the 
complex domain; failing to notice this may cause an erroneous conclusion. 

Remark For a fixed z, the presence of the functions A and A in a/ and rji 
themselves only affects their signs. More specifically, the only dependence on A and 
A takes the following forms: 



1 AHQ 1 AHS) 

F{z) = -F2" + {F^f + 3F2F2' - 2F|F2 + 2C'^F^ + 1^1(3^2' - 4^2 F2). 

Hence the product is a function of z only. According to Cartan [15], this sign 
situation is accounted for by a local CR diffeomorphism and therefore does not generate 
a new CR structure. Hence we have proved the following theorem. 

Theorem 3 Locally, the CR structure (111) (as an equivalence class) determined by 
the function c given in fll8|19p is independent of the choice of the function A{() 7^ 0, 
once the form of F2{z) is fixed. 

Altogether, the freedom of choosing various A{Q 7^ does not affect the CR 
structure of a type N metric which we are considering. Hence, for the simplicity of 
representing new metrics distinguished by CR structure, we can just set A{() = A{Q = 2 
(see the Conclusions). In hindsight, the classical symmetries we have obtained are 
nothing more than a particular manifestation of the underlying CR equivalency. We 
believe this connection between the two may as well suggest a more general concern if 
one aims to find, through the (classical or higher) symmetries, additional exact solutions 
to the Einstein equations formulated with CR structures. 
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We will see later examples of solutions that have constant CR invariants, and 
remarkably, one of them is the solution of Leroy-Nurowski. Nonetheless, this feature is 
generally not true for other solutions. 

7. Conformally flat solutions 

Before we try to solve fl2T|) for type N solutions, it is important to find out in advance 
those conformally flat solutions satisfying \E'4 = which are automatically contained in 
the general solution of ( 12T1) . We insert the ansatz fll8lll9p into the expression for \l/4 
given by ( !T7|) . and re-normalize to pull out just a simple complex- valued function of 
z: 

4Acos-^ 

= 2F,F[' + Q{F[)^ - 10(F2 + iCi)F,F[ + {-F^ + 3F| + 61^1^2 - 3Cl)Fl (29) 
We now apply ( l20l) and use (12T]) to substitute for J'", which gives us 

K = [Ajr - |A(J')' + 2(AV2 - 2C^)f] + i [-2Ci {J" + 3AJJ')] , (30) 
or in terms of P{J), 

K = P [UP' - lAP + 2AV2 - 4Cl] + i [-2C1P (P' + 3A J)] , 

where we have put the real and imaginary parts in separate brackets. Replacing J" with 
the help of ( 130!) . we can rewrite the equation (12T!) as 

= lAK{J'f - {2AKJ + K'){AJ - 2iCi)J' + 

which clearly has K = 0, i.e, all conformally fiat solutions, as some of its solutions. 

If P{J) is not restricted to the real domain, then solving the first-order ODE K = 
for P(J) leads to the following general solution 

P(J) = -|a(j^ + M)+c.(j±H^)"'. (31) 

with a complex constant C2. 

If, instead, we restrict P(J) to be real, a simultaneous vanishing of the real and 
imaginary parts of K respectively yields the following set of two equations, provided 
P^O, 

= Ci(P' + 3AJ), 

P' = ^-2AJ+1^ 
3 J + A J ' 

both of which are consistent with fl22l) . There are now two cases for solutions. 

The case Ci ^ requires that both ODEs be satisfied, so that we have a unique 
solution 

P(J) = -hj^ - ^ (32) 
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which, by solving J' = P{J), gives rise to 

Atan(3Ci(^ + Co))' ^ ' 

In the hmit Ci 0, the above solution becomes even simpleiEl: 

'=^A(lW- '''' 
From (l33i) . we have 

F = + v^^i F = 

' ssin(3Ci(2 + Co))' ' tan(3Ci(^ + Co)) ^ ' 

with negative- valued A = — s^. Note that it is only at this stage that the reality 
condition on Fi 2, i.e., J' > 0, requires A < 0, i.e, a negative cosmological constant. An 
important remark that can be made is that the extended form of the Leroy-Nurowski 
solution (see the next section) resembles this solution greatly with simply differences in 
the coefficients. 

For the other case when Ci = 0, we have 

P(J) = -|AJ2 + C2J^/\ (36) 

with a real constant C2. From f l23|) . the solution J{z) is determined by 

1 



dJ = z + Co. (37) 



2 

Since J' > 0, we cannot have both A > and C2 < 0. Hence, we can discuss three other 
sign possibilities, the details of which are put in Appendix C. 

We note that the special solution (IMI) corresponding to Ci = C2 = serves as the 
single "point" where these two families of conformally flat solutions are joined up. 

Modulo possible sign differences in aj and rji caused by square roots as already 
discussed, the Cartan invariants for both (133|) and (!34|) . as calculated via (l28l) and the 
equations for the other invariants, as presented in Appendix B, are given by 

414/2 . 41 29 



e \ 5 2VT0 2VT0 ^^^^ 

, /2 i 219/4 327 

Ut = 6i\ -, riT = — ■ , Ci = , e = ±1. 

Remarkably, they are all constant and do not depend on Ci. Nonetheless, this is not the 
case for the other conformally fiat solutions obtained from f l57|l with C2 7^ of which 
the Cartan invariants are generally functions of z and C2. For instance, simplified by 
( 136|) and J' = P{J), the first Cartan invariant satisfies 

2. /2 /3Aj4/3 + 2C2\^ , , 

a;(..C,) = -16^-(3^jj^j , (39) 

where J = J{z) belongs to one of the three cases described in Appendix C. 



^ Both (|33|) and (p4| would be of particular importance for the perturbation theory on type N 
solutions near flat ones. 
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Two conformally flat Einstein spaces may have non-equivalent CR structures. This 
does not conflict with the previous theorem because in a conformally flat spacetime, one 
is free to make different choices from among the multiple shearfree null congruences and 
therefore may have non-equivalent CR structures attached to them. 



8. An extended form of the Leroy-Nurowski solution 

Now we can reveal a fuller extent of the exact twisting type N solution first discovered 
by Leroy, and re-derived by Nurowski within the framework of CR geometry, upon the 
latter of which our current work is mainly based. We hope that our derivation of this 
solution will make the process behind the previous discoveries appear clearer. 

Given Nurowski's form of the solution (see [11] or ( H7|) ) and recasting it into the 
form of the ansatz (118|ll9p and ( l20l) . we find the following special solution to ( l22|) 



P(j) = _1aj2_M (40) 
^ ^ 3 4A ^ ^ 

which gives rise to a solution to fl2Tl) : 

J = ^ (41) 

2Atan(|Ci(^ + Co))' ^ ' 

In the limit Ci — ?■ 0, the above expression becomes even simpler: 
A{z + Co) 

which is quite similar to that of f l5^ . Back to the case with Ci ^ 0, using (^U^, we have 

' 2ssin(iCi(^ + Co))' ' tan(iCi(^ + Co)) ^ ' 

with a negative A = — s^. Note that it is only at this stage that the reality condition 
on Fi requires A < 0. In the end, our extended version of the Leroy-Nurowski solution 
takes the form 

P(C, = ± (44) 

2ssmh{^Ci{z + Co))y/AA 



2Ci 

tanh(iCi(2 + Co)) 



+ ; . ^ + Ci 



(45) 



i.dc - [ idc] . (46) 



A ^ J A 

The flexibility of choosing the function A{Q and real constant Co,i may perhaps facilitate 
a possible future application of the solution. From this extended version, one can obtain 
the original form of Nurowski [11] by setting 

AiC) = CiC, AiC) = CiC, Co = 0, 
and consequently. 
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Note that all Ci's are canceled out in the above expressions. Hence another way of 
obtaining (HT]) is by taking the limit Ci in (jH]) and (jl5]) (cf. (jl2])) and setting 
Co = and A{C) = 2. 

Modulo possible sign differences in «/ and rji caused by square roots as already 
discussed, the Cartan invariants calculated from (jUj) and (H2l) are both given by 




1 2^/2 I 
^/ = i\/7, ^^ = -i - 31/4 -53/4 ' ^^ = -^' " = ±^- 
Like (138!) . they are all constant and do not depend on Ci. 



(48) 



9. An example of power series solutions 

For simplicity, assume that Ci = in the ODE fl2T]) . Now consider the power series 
solution of ( I2T]) satisfying the regular initial conditions J(0) = 0, J'(0) = uq > and 
J"(0) = 0. A simple calculation gives us the first few terms of this series 

00 

J{z) = J2 ^i^'^^ = ^oz - lAulz^ + ^A^ulz^ + ■■■■ (49) 

1=0 

Moreover, this series solution, convergent in a neighborhood of 2 = according to the 
Cauchy existence and uniqueness theorem, is of type N with a non- vanishing Weyl scalar 
oc K{z). Particularly, 

To see that the solution fH9|) is not equivalent to the Leroy-Nurowski solution, we 
calculate the first Cartan invariant ai{z) via (128|) which, in this case, is no longer a 
constant. In particular, this series solution has 

0/(0) = 0, 

with ai{z) and also K{z) continuous at z = 0, while the values given in (l48ll are always 
non- zero constants. This is sufficient to assert that the ODE (1211) as well as its reductions 
fl22|) and fl26|l indeed contain new twisting type N solutions. 



10. One-parameter deformation from a conformally fiat solution to the 
Leroy-Nurowski solution 

One feature that makes fl22|) preferable to the other two ODEs fl2Tl) and (|26|) is that 
the conformally flat solution fl32|) and the extended Leroy-Nurowski solution pUj) are 
just simple quadratic functions, without poles in the complex plane, compared to their 
counterparts (133|) and (14T|) . Also note that these quadratic solutions with Ci = 
do not correspond to any solution of the Abel equation (l26l) since the form (1271) with 
the non-constant function 2'(J), excludes all quadratic functions as solutions. These 
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well-behaved quadratic solutions facilitate a study of the power series solutions near 
them. Additionally in Appendix E, we also present a Puiseux series solution to f l22|) . 
the existence of which is suggested by the weak Painleve tests performed in Appendix 
D. 

To simplify the notation, we apply the scaling transformation^ J = Ciw/A, 
P(J) = Cfg{w)/A with A 7^ 0, Ci 7^ such that ([22]) takes on the form already 
noted as ([1]) with C = 1: 

^ _to- + 2uf _ 2 _ 10 

29 9 3 ^ ' 

We look for power series solutions for this equation corresponding to the regular initial 
conditions g{0) = Mq 7^ 0, g'{0) = 0. The first few terms of this series read 



j=0 



5 (no + I) 2 2 (no + |) + 6) ^ 76 (no + |) + 6) (no + 1)6, 
"° 27u^ " 1215^ " + 



(51) 



where all odd order terms vanish. The remainder of the coefficients in the series can be 
determined by a recursion relation which is valid beginning with uq. 

k-l 

= {2k + l){k + l)uoU2k+2 + {2k + I) U2k + J2^k + /+!)(/ + l)u2i+2U2k-2h k>2, (52) 

1=0 

while U2 and ^4 can be easily read off from (1511) . It is clear that this relation allows 
one to calculate the coefficients to whatever order desired. One can easily see that the 
coefficient of w'^'', namely U2k, is a /cth-order polynomial, Pk{uo), divided by ul^~^ 7^ 0. 
Remarkably, this infinite series reduces to simple quadratic functions in two special 
cases. The reason for this is that for every value of /c > 2, the polynomial Pk{uo) has 
the factors (no + |) (no + 6), as can be seen in the few terms demonstrated in (1511) 
above and can easily be shown by induction. Hence for no = — |, we retrieve the 
Leroy-Nurowski solution (HOl) . which in this notation is simply 

9LN = - ilw^ + f ) . (53) 

As well, for no = —6, we retrieve a conformally flat solution fl32|) . which has the form 

9cF = -{lw' + Q). (54) 

For all other values of no 7^ 0, the formal series solution ( ISTj) may then be viewed 
as a generalization of these two known solutions, in terms of a power series with 
inflnitely many terms. It is interesting that in every one of these polynomials, Pk{uo), 
all coefficients are negative, so that the only possible real roots would be negative. Our 
numerical calculations suggest that none are smaller than —6, and that there are no 
other roots common to all these different polynomials. 



Once having a solution g{w), one may ehoose a sign for A in order to have P{J) > 0. 
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The series (I^Tj) does define, in the complex domain, a function holomorphic in 
some neighborhood of the origin as is shown by the following method of determining a 
non-zero radius of convergence for it. We present the proof in Appendix F. 

Theorem 4 Given the series ( 15T]) with the recursion relation ( 152|) and a fixed 
no 7^ 0, one has the following bound: 

|^2il < -^2^, J = 2,3,---, (55) 

provided that one can pick two constants C > and M > such that they satisfy 



2 (mo + f ) (^^0 + 6) 



Tlul 



16 

5 1 \ 9 /tt^ 1 



< ^. (56) 



■ , , / o ■ . ^C < luol. (57) 

_3 \uq\) 4M2 Vl2 4y - ' ' ^ ' 

The existence of such an upper bound fl55l) on guarantees a lower bound M^^ 
on the radius of convergence. For instance, if we take uq = —2, which lies nicely in the 
interval between — | and —6, we can at least pick 

satisfying both ( l56l) and ( 1571) . The bound (l55ll is by no means optimal at every uq ^ 0. 
In fact, our numerical integrations of fjSOj) with uq sampled between —6 and — | all 
indicate that in the real domain, the series solutions ( ISTl) with —6 < mq < — | are all 
well sandwiched between the parabolic curves of flS^ and flM|) . and therefore suggest an 
infinite radius of convergence on the real line. Moreover, by applying the transformation 
w — )■ ^ to (|5Ui) and studying the formal (Puiseux) series expansion of the transformed 
ODE at the origin, we find the following asymptotic expansior0 of (|50|1 as — )■ oo (cf. 

my 

g ~ -\w^ - 6 + Ui/^w^'^ + 0(u;^^/3), 

where M4/3 is an arbitrary constant. This asymptotic behaviour at infinity, consistent 
with our numerical calculations, again suggests that we may significantly extend the 
radius of convergence for flSTl) at least in the real domain. 

An additional comment is that the Cartan invariant a/, computed from ( ISTl) with 
—6 < Mo < — |, is not constant and has a dependence on Ci, contrary to the special 
cases for those values of uq at the two endpoints of the interval of values for uq being 
considered. 

11. Conclusions 

We have begun with the advantage of prior work done on the use of (3-dimensional) 
CR manifolds to look for solutions of the Einstein field equations that correspond 



^ We also find another asymptotic expansion that has the first two leading terms identical to (|53p . 
but also involves fractional powers of w in a complicated way, hence not presented here. 
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to algebraically-degenerate Einstein spaces, with twisting principal null directions. A 
general solution of those reduced field equations for the two functions of three variables 
would generate all twisting solutions of Petrov type N. Of course we did not achieve 
this; however, after the assumption of a single Killing vector in a particular direction, 
our ansatz for group-invariant solutions obtained from the infinite-dimensional classical 
symmetries of the field equations, allowed us to obtain a single ODE, the solutions 
of which would generate a family of solutions of type N with twisting principal null 
directions. That ODE is either a rather simple, third-order nonlinear equation for 
J = J{z) in which the independent variable z does not appear or, equivalently, an 
even simpler, second-order nonlinear equation for g = g{w), where w is a. dimensionless 
re-scaling of J and g is a. re-scaling of J', which includes a non-zero value for A, the 
cosmological constant. Within the same ansatz, we have also investigated all the cases 
of solutions corresponding to conformally fiat spacetimes to which type N solutions may 
degenerate, which helps us look for non-trivial cases. 

We have studied this second-order equation at some length. In particular it contains 
one parameter, Ci, which may always be re-scaled to the value unless it happens 
to be zero. In the case that it is zero, the equation can be reduced still further, to a 
first-order equation of Abel type. Following standard approaches to Abel equations we 
were unable to determine any method that we thought would generate reasonable type 
N solutions, although this is still an ongoing project of considerable interest. However, 
when Ci is not zero we have considered various sorts of solutions which it might have. 
We have shown that it does have solutions which are holomorphic, in the complex plane 
in a neighborhood of the origin, and have found an asymptotic behavior near the (real) 
infinity. In particular we have picked out especially those solutions which are even 
functions of w and looked at power-series solutions about the origin, both analytically 
and numerically via a Maple computer program. We have determined a moderately- 
simple recursion relation for the coefficients of the powers of w"^ in the series solutions, 
which determines the coefficient U2k+2, of w^^"*"^ {k > 2) in terms of all the previous 
coefficients, looking at all of them as determined by the value of g{0) = uq. This series 
terminates quickly for just two particular values of Uq, in the form —a{uo)w'^ + Uq, 
with a constant, different for the two values of uq. The value mq = — | generates the 
previously-known Leroy-Nurowski solution, while the other one uq = —6 is unfortunately 
simply a conformally fiat solution. To ensure that these series solutions are distinct 
from the Leroy-Nurowski solution, we have used the work of Cartan on the question 
of the equivalence of two CR manifolds, which requires the equality of the set of six 
Cartan invariants. We have found that any value of uq between these two special values 
generates Cartan invariants that are quite different from those at the endpoints of this 
interval, and therefore distinct from those of the Leroy-Nurowski solution. 

The solutions characterized by values of mq between —6 and — | have an asymptotic 
behavior, via a Puiseux series around the (real) infinity, that has the same form 
— a(— 6)w^ — 6 as the conformally fiat solution aforementioned, but also lower-order 
terms involving third-roots of w, which doubtless generate algebraic singularities there. 
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Numerical integrations via Maple agree with this behavior, showing negative values of 
g{w) as needed and very simple structure, for all real values of w. The same numerical 
integrations do show singularities in the solutions for mq > — |. As well, numerical 
calculations of the coefficients U2k+2, for several values of Mq E (—6, — |) (e.g., uq = — f^^) 
show that starting at a large enough k, they alternate in sign while their absolute values 
are monotonically decreasing at rapid rates. We therefore postulate that these solutions 
are everywhere non-singular and well-behaved on the real axis, and believe that they 
might define new well-behaved, transcendental functions with algebraic singularities off 
the real w-axis. The proof of such a conjecture is still being pursued; nonetheless, we 
feel that the numerical calculations justify the belief that this is a sufficiently interesting 
result as to merit the attention of a wider audience. 

To conclude the discussion, we present here our new class of metrics, which, without 
loss of generality, may be considered by setting A{() = 2 in the ansatz ( ITSi) . Although 
we will present it here with the new real coordinate z introduced in ( IT9|) . with this choice 
of A{() it is the same as the usual coordinate y used in (I9l). As well, our studies with 
the equation for P{J), equivalently g{w), allow us to replace z by its form in terms of 
J as determining the imaginary part of d^, via d^ = dJ/P{J), namely, 

( = x + iz = x + iz{J), d( = dx + idz = dx + —dJ. 

For simplicity of presentation, we show both forms below, with coordinates {x,2;,u,r} 
or {x, J, M, r}: 

g = o [dCdC + A (dr + WdC + WdC + H\)] 

COS 2 / 

with real- valued J = J{z), J' = dJ/dz = P{J) > and P' = dP/dJ such that 

^=1{^ + ^J + i^i) (^"" + 1) = i ilP' + ^'^ + i^'"' + 1)' 

H = -|A/cos(r) = -|APcos(r). 

where Ci is an arbitrary real parameter. The function L as in d = — Ldu can be 
chosen so as to be real-valued: 

L = -e"^^" / exp ( / F2dz] dz = -e'^i" ! T^^^vi I ^ dJ] dJ, 



such that from (Q, 

eCi^du - 2 [/exp (/ Fsdz) d^] da; _ e^^^du -2[J p-^ exp (/ FgP-Mj) dJ] dx 
~ exp (/ F2dz) ~ exp (/ Fa^-^dj) ' 

where F2 is given by 

F2 = ^-AJ=lP'-AJ. 
Meanwhile, the functions J{z) and P{J) respectively satisfy 



( T")"^ in 

r' = S±L _ 2AJJ" - yA(J')' - 2(AV2 + J', 

^„ ^ jP' + 2AJ)2 _ ^ _ 

2P P 3 ■ 
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In particular, the original metric by Nurowski [TT] corresponds to the case Ci = 0, 
J = and a proper choice of the integration constants in A. 

We can here note the philosophy that certain ODEs themselves may serve the 
purpose of defining new transcendental functions; for instance, we recall the Painleve 
functions and the associated ODEs. Hence our situation with new type N solutions being 
determined by a second-order nonlinear ODE is presumably not too different from that 
of Mauser's solution (in terms of hypergeometric functions [1]) which is determined by 
a second-order linear ODE, although it is true that there has already been much more 
extensive studies made on the properties of hypergeometric functions than have been 
made for newer functions defined by solutions of nonlinear ODEs that may not even 
have the Painleve property. 

Appendix A. The Abel ODE 

The equation (l26l) actually does have the following special solution 

^- = -¥T6- 

However, it can be shown to correspond to a conformally flat solution ( ]36|) . and hence 
is not interesting. 

Unfortunately, we have had no luck so far finding the general solution to ( 126|) or 
any other special solution other than (lA.ip . Since constructing the general solution 
to the generic Abel ODE has remained an open problem for decades, the general 
strategy of integration nowadays mainly lies in recognizing, within a suitable class of 
transformations, the ODE in question as equivalent to a previously solved equation. 
Such a procedure has been programmed into the current state-of-the-art Maple code 
dsolve (or abelsol) [T8l[T9], which presumably covers all/most of the integrable classes 
presented in Kamke's book [26] and various other references (e.g., [27]). However, this 
code, as tested by us, does not recognize ( l26l) as a known solved type, e.g., the AIR 
class. Other attempts by us, such as the symmetry method, on finding special solutions 
all have failed or just led to (lA.ll) . 

So far we have not been able to find a similar reduction for the ODE fl22|) with 
Ci 7^ 0, nor can we negate the possibility that fl22|) with Ci ^ may contain different 
type N solutions other than the case with Ci = 0. In fact, the Cartan invariants 
calculated with (12T]) generally do have a dependence on the constant Ci even though 
this is not the case for all the conformally fiat solutions and the Leroy-Nurowski solution 
(see dSHD and (gHD). 

Appendix B. Cartan invariants 

A concise description of Cartan invariants can be found in Section 2.1.3 of [21]. Given 

c = c(C, C) and 

r = I (<9^^+ 2cl) , / = —di^d^c — cd^c 
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where aj is presented, the next three Cartan invariants, when r 7^ 0, 

1 



taken from 
read 



/3/(C,C) 



32 (rf) 9/4 



3r d^r (9^r + 3r d^^r d^r — rr i dc^r d^r 



7/(C,C) 



7d(^r d^r + lQcrd(^r + IQcrd^^r — Srrd^c + IGccrr 



Tr'^d^r d^r + Ir'^d^r d(^r — rr ( 8rd(^d^r + 8rd(;d^r 



32 (rf) L 

+ d^f d^r + d(^r d^f + Acfd^r + Acrd/^f + Acrd^f + 4cf(9^r 
+ 24rf9^c + 16ccrf j 



^/(C,C) 



16r(rr)'^/4 



— 2d^rd^r — Acrd^r — Acrd^r + IQrrd^c 



Due to the length of C,i and r// as calculated with Maple for our studies, we will not put 
them here. All Cartan invariants are uniquely determined by the function c = c(C, C)- 



Appendix C. Conformally flat solutions 



For a further integration of fl37l) . we have three separate cases. 

Case i; A < 0, C2 > 0. We always have J' > 0. Then the solution is determined 

by 

In ^ + 2 arctan — = -2V2M^z + Cq , 

G^-V2GM + M^ \M^-G^ ^ 



M = ± 



2C2AV3\1/^ 



G = (AJ) 



1/3 



In the real domain, the inverse function J = J{z) is well defined over z + Co G 



V2|A/3|' V2|Af3|) 



instead of the entire real line, and has singularities at z + Gq 



-V2\A'P\' 



Case 2: A < 0, C2 < 0. We need |J| > (2C72/3A)^/^ for J' > 0. The solution is 
determined by 



In 



M + 


G 






G 




M - 


+ 2 arctan ( ^ j 



2M^(z + Co), 



2G,^"' 



3 J ' V . ' ' ' - V 3A 

In the real domain, the inverse function J = J{z) is well defined over z + Gq & 



-00, 



■2|A'/3| 



U 



(21^^' ' singularities at z + Cq 



-2\M'^\- 
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Case 3: A > 0, C2 > 0. We need |J| < (2C2/3A)^/^ for J' > 0. The solution is 
determined by 



M + 


G 


fG\ 


M- 


G 


+ 2 arctan — 


\m) 



= 2M\z + Go), 



In the real domain, the inverse function J = J{z) from above is well defined over the 
entire real line. 

Appendix D. Weak Painleve tests 

We have three ODEs f l2T|) . f l22|) . and f l26|) at hand that may be explored for new 
twisting type N solutions. A particular, probably useful way to decide which one of 
these equations has a better chance for one to find a solution is given by the (weak) 
Painleve test [291 EOl Elj. This test reveals the nature of the movable singularities 
of the general solution of a nonlinear ODE. Failing the test means the occurrence of 
certain undesirable movable singularities, e.g., infinitely branched singularities, that 
relate to non-integrability [31], although it may still be possible to find special solutions. 
Associated to the (weak) Painleve test is the global property called the (weak) Painleve 
property. An ODE possesses the Painleve property if the general solution can be made 
single-valued. For the weak Painleve property, it requires that the general solution be 
at most finitely branched around any movable singularity. The tests themselves are by 
design sets of necessary conditions respectively for these properties. 

In this appendix, we will show that none of the three ODEs pass the Painleve test, 
and that ( |2T1) also fails the weak Painleve test while the other two pass. To begin, we 
detail the test procedures on fl22|) . Then we briefly comment on f l26|) and simply point 
out where the tests fail for f l2T]) without dwelling on details. 

The equation (l22l) surely does not have the Painleve property for the coefficient of 
the {P'Y term clearly violates the necessary conditions for the Painleve property [29] 
(see p. 127). This is also confirmed by the test conclusion that ( l22l) has movable algebraic 
singularities. 

Step 1 [Dominant behaviours). Assume the leading behaviour of a solution P{J) 
to be 

P ~ uox"", x = J - Jo, Uo^O, m^O, 

with m not a positive integer. Substitute this form into fl22|) and select out all possible 
lowest order terms as listed below 

Since m 7^ 1, we only have two possibilities. For m < 1, x^*""^ is the lowest order term 
and the vanishing of its coefficient requires 
2 
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given UQ,m ^ 0. For m > 1, the constant 2(A^Jq + Cf) is the lowest order term, which 
does not vanish in general, hence not interesting for the purpose. To summarize, we 
obtain m = ^ with arbitrary uq ^ 0, i.e, that 

is the only detected dominant behaviour. 

Step 2 [Resonance conditions [29] (see p. 87)). Having found the dominant 
behaviour, now we consider the possibility to extend it to a Puiseux series expansion 



,^ja+2)/3_ 



p = ^n,(J- Jo 

j=0 

This requires the determination of the locations (j + 2)/3, called Fuchs indices or 
resonances, where arbitrary coefficients may enter the Puiseux series. Consider the 
dominant terms 



E{J,P)=PP" + \{P'f 

of that contribute to the leading behaviour x^"^~^ = X^^'^- Then compute the 
derivative 

The Fuchs indices satisfy the so-called indicial equation 

Hence we obtain a fractional resonance at (j + 2)/3 = I with J = 0. 

Step 3 (Compatibility conditions). At j = 0, we know, from the the first step, that 
Uo{^ 0) is indeed an arbitrary coefficient. This completes the test. In conclusion, ( 122|) 
passes the weak Painleve test. 

Remark Note that no pole is detected from the test above. The ODE for P^ still 
involves a Puisuex series instead of a Laurent series since the cubing does not eliminate 
all third roots of x- According to pT], The presence of movable algebraic singularities 
is not incompatible with integrability. 

The very design of the weak Painleve test limits its usage only as necessary 
conditions for the weak Painleve property. The test can neither detect movable 
(branched) essential singularities themselves nor exclude an accumulation of algebraic 
singularities forming a movable essential one that may be severely branched. These 
possibilities make a rigorous proof of the weak Painleve property not at all a trivial one, 
which by itself may deserve a specialized article to discuss. See examples in [321 [331 iM] • 

According to Painleve [351 [36] , the only movable singularities of solutions to the 
first-order ODE y' = F{x, y) where F is rational in y with coefficients that are algebraic 
functions of x, are poles and/or algebraic branch points. In addition, the only nonlinear 
ODE in this class that has the Painleve property is the Riccati equation which ( 126|) is 
certainly not. Hence the equation (126|) automatically has the weak Painleve property, 
but not the Painleve property, and it is free from movable essential singularities. 
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The equation (1211) admits two families of dominant behaviours (cf. fl34p and 

2 4 7 

J ~ —— r, Fuchs indices = —1, -, -: 

3A(2-^o) 3' 3' 

J ~ -TT— Fuchs indices = —1, 

A(^-^o) ' 2 ' 2 

It fails the weak Painleve test for having irrational resonances. This means that fl2T]) 
has an infinitely branched movable singularity, which is a strong indicator for non- 
integrability [3Tj . 

Since our attempt of solving (l26l) has not been successful, we decided to focus on 
f l22|) and explore some of its features that may facilitate constructing new solutions. 



Appendix E. Puiseux series solutions 



As indicated by the weak Painleve test, the ODE ( 122|) for P{J) possesses a formal 
Puiseux series solution 



k=0 

= no(J- Jo)'/'-3AJo(J 



9(A2J2 



h 

2 



20no 

27(109 A2J2 + 36C2) (A 



2 72 





2800^3 



(J 



(E 



with two arbitrary complex constants Mq 7^ and Jq. In particular, this Puiseux series 
solution contains a special case for Jq = ±2iCi/A (A 7^ 0) such that 



P = Uo{J± 



2iCi 



2/3 



^A J2 + 



4C| 
A2 



This finite expression coincides with the known solution ( pTl) (setting uq = C2). 

Theorem 5 Given that mq 7^ and uq, Jo G C, the ODE f l22|) admits a formal 
Puiseux series solution f lE.ip such that it converges in a neighborhood of Jq. 

Proof: The idea of the proof, following many standard proofs of Painleve property, 
is to convert the Puiseux series into a power series solution of a regular initial value 
problem (e.g., [321 El])- First we define 

Z = P^/^{P' + AAJ). (E.2) 

Then differentiate it once with respect to J and substitute P" using fl22l) . Hence we 
obtain 

2(AP -3A^J^ -3C?) , 
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The system (1E.2|E.3P is equivalent to the ODE 
U = P^/^, we can transform the system into 

dJ _ 2U^ 

dU ~ Z-AAUJ' 

dZ _ 4(3AV2 - At/2 + 3Cf)U 

dU ~ 3{Z - 4AUJ) ■ 

which has a unique power series solution about U 

2 



. Now by introducing a new variable 



(E.4) 



J=Jo + 



3^0 
2(A2J2 



CI 



(E.5) 



By the Cauchy existence and uniqueness theorem, both series have non- vanishing radii 
of convergence. From the series (lE.Sp . the corresponding solutions to flE.2|E.3P then 
take the form 

2/3 



3^n 



{J - Jo) 



k=l 



z = Zo + J2MJ-JoY'-''^^'- 



k=0 

with Zq ^ 0. This completes the proof. 

The series (lE.ip clearly contains type N solutions that are not equivalent to Leroy- 
Nurowski's since they all continuously deform to the conformally flat solution (15^ in 
the limit Jo — )■ 0, Ci — )■ 0. We already know that the latter has a non-constant Cartan 
invariant aj given by f l5^ . 



Appendix F. Proof of Theorem 4 



The induction begins with 



F4I < 



16 



which holds by the assumption ( !56l) . Now assume that for > 2 and j = 2, ■ ■ ■ , k, the 
bound fISS]) is true. Then for A; > 3 and 1 < / < A; — 2, we can bound the product 

U2i+2U2k^2i by 



\U2l+2U2k-2l\ < 



< 2 
= 2 



(2/ + 2)2(2A;-2/)2 
(2A; - 2/)2 + (2/ + 2)2 



(2A; + 2)2 
1 1 



(j2^2k+2 



(2/ + 2)2 {2k-2iy 



(2/ + 2)2(2A;-2/)2 

C2]Vf2fe+2 



{2k + 2 



\2- 



(F.l) 



(2fc+f) 


U2k 


+ (2P + A:+ l)|M2M2fc 


, Y!l=l{k + I + + l)\u2l+2U2k-2l\ 


(2A; + l)(A; + l)|no 


(2A; + + l)|uo 
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The second inequality above is due to (a^ + 6^)/(a + 6)^ > |. Rearranging (1321) and using 
the triangular inequality together with fIF.ip . we obtain an upper bound for |u2fc+2h 

\u2k+2\ < 

^ {2k + I) + {2P + k + 1)\U2\ CM^'' S{k) C^M^''+^ 

- i2k + l){k+l)\uo\ (2A;)2 ^ 2{2k+l){k + l)\uo\' i2k + 2y' - ^' 

(F.2) 

where we define 

We can evaluate the first summation above in terms of the digamma function 



4-4 (/;; + / + !)(/ + 1) ,^ /,T/,N 



where 7 is Euler's constant, which is approximately 0.57721 • ■ -. The second summation 
has the following bound 

'^{k + l + + 1) _ (2A; - l){k - I) 



{k-ly ^ + 

<2.^g^^(^-2).^-2P.(i,.)<(^_2).^. 

Note that the trigamma function \E'(l,/i;) > for all integers k > 3 and that 
\E'(1, k) ~ A;~^ for k — +00. Combining these two bounds, for k > 3, we obtain 

S{k) ^ (7rV3 - 2) P + k'^ik) 



2{2k + + 1) ~ 2(2fc + 1)(A; + 1) 

(7rV3 - 2) P + A;2 yr^ 1 



4P 12 4' 

where we use the fact that < "^{k) < k for all integers k > 3. In addition, the first 
term in (\F.2\\ is bounded by 

(2fc + |) + (2P + fc + l)|n2| CM'' ^ + I) + + ^ + 1) (1 + m) CM^^ 



(2A; + 1)(A; + 1)|mo| (SA;)^ ~ (2A; + 1)(A; + 1)|mo| {2k)' 

{k + 1)2 CM"'+' 



< 



< 



\ 


8< 


1 






\uo 


\ 




1 









k'M' {2k + 2)2 
9 CM2'=+2 



4M2 (2A; + 2)2' 

where the last inequality becomes an equality for k = 2. Altogether, we obtain for k >2 



1 

M2fc+2 < 1 — r 

Pol 



6 1^5 +S'M-\^C 



3 \uo\J 4M2 V 12 4 



CM2'=+2 CM2'=+2 
< 



(2A; + 2)2 - (2A; + 2)2' 



given the assumption f l57p . Here Si is the Kronecker delta. This completes the induction. 
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